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Resumen 


Tanto la agricultura como las descargas de aguas residuales son fuentes importantes 
de nitrógeno en las corrientes de agua superficiales. El objetivo en este trabajo es 
describir de forma espacial y temporal las concentraciones de especies de nitrógeno 
en el río Zahuapan y sus relaciones con algunos atributos de las microcuencas, como 
área, pendiente, cubierta vegetal, agricultura, erosión, población y precipitación 
pluvial. Con el programa TAS se definieron los polígonos de las microcuencas del 
río. Los resultados muestran que las concentraciones de las especies de nitrógeno 
disminuyen en las microcuencas que se encuentran en las cercanías de la presa de 
Atlangatepec. De manera temporal, las concentraciones tuvieron su mayor valor en 
los años de 2001 y 2002. Las concentraciones de NT, NTK, N-NH, y ss aumentan 
en los meses de estiaje. Por otro lado, las concentraciones de N-NO, y N-NO, 
disminuyeron en esta temporada del año. La relación entre las concentraciones y la 
precipitación mensual resultó tener el mayor número de coeficientes de correlación 
significativos que los realizados con la precipitación diaria. Los coeficientes de 
correlación entre las especies de nitrógeno fueron significativos a un nivel de 0.05. 
El análisis de regresión múltiple de la forma N-NO,= A, + A,*NT + A,*N-NH, + 
AN ,, mostró coeficientes mayores que 0.5 en nueve de las doce microcuencas. Se 
concluye que la presa de Atlangatepec juega un papel en las concentraciones del 
nitrógeno en el río Zahuapan. 
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Introducción 


El ser humano ha impactado el ciclo del ni- 
trógeno debido a que el exceso de fertilizante 
nitrogenado que se emplea en la agricultura 
es transportado a los cuerpos de agua o se 
infiltra a las aguas subterráneas. Además, las 
aguas residuales sin tratamiento contienen 
nitrógeno y son descargadas a ríos o embalses. 
El agua residual municipal cruda puede tener 
concentraciones promedio de nitrógeno total 


de 35 mgL” (Hammer y Hammer, 2001). El 
río Zahuapan (Tlaxcala, México) ha sido el 
receptor de los escurrimientos y las descargas 
de agua residual cruda o tratada de los 
centros de población de la subcuenca, por 
lo que resulta importante saber cómo han 
cambiado con el tiempo las concentraciones 
de las formas químicas del nitrógeno en el río 
Zahuapan y cómo tales concentraciones están 
relacionadas con las fuentes. La relación entre 
las concentraciones de nitrógeno en los cuerpos 


Tecnología y Ciencias del Agua, vol. IV, núm. 4, septiembre-octubre de 2013, pp. 47-61 


Ciencias del Agua, vol. IV, núm. 4, septiembre-octubre de 2013 


« vj30]0u32], 


Muñoz-Nava ef al., Relación entre el nitrógeno y las microcuencas del río Zahuapan, México 


de agua y las fuentes varían de acuerdo con 
las condiciones ambientales, pero las rutas 
del ciclo del nitrógeno son las mismas. Por 
eso los estudios realizados en otras latitudes 
son de gran utilidad, pues muestran las 
relaciones entre estas variables. Por ejemplo, 
Hayakawa et al. (2006) encontraron una 
correlación positiva entre las concentraciones 
de nitrato en una corriente de agua al este de 
Japón y el área que se encuentra en la parte 
alta de la cuenca de drenaje, aumentando 
las concentraciones en cuencas con uso de 
suelo urbano, comparado con aquellas donde 
hay humedales naturales y bosques. Una 
relación específica la proporcionan Tesoreiro 
et al. (2009), quienes señalan que el agua 
subterránea es la fuente principal de nitrato 
en cinco corrientes superficiales de Estados 
Unidos. Puckett y Hughes (2005) reportan 
una relación entre las concentraciones de 
nitrógeno y las precipitaciones pluviales, 
al encontrar que las concentraciones de 
nitratos y amonio, en un terreno de cultivo 
sobre un acuífero superficial en Carolina del 
Sur, Estados Unidos, aumentan en periodos 
de mayor flujo de la corriente. En cuanto a 
una estimación cuantitativa del aporte de 
nitrógeno, los cultivos de maíz y soya son los 
mayores contribuyentes, con 52%, seguido 
de las deposiciones atmosféricas, con 16%, 
en las cuencas del Mississippi y Atchafalaya 
(Alexander et al., 2008). En Inglaterra y Gales, 
el aporte de la agricultura es de 60% y de las 
descargas de aguas residuales de 32% (Burt et 
al., 2011). En particular para México, se han 
realizado estudios que relacionan el nitrógeno 
con fuentes puntuales y difusas. Izurieta et al. 
(2002), empleando un modelo de funciones de 
carga, estimaron en la cuenca del lago Cuitzeo, 
un porcentaje de contribución de nitrógeno por 
fuentes puntuales y difusas de 66.6 y 33.4%, 
respectivamente. En un estudio realizado por 
Silva et al. (2002) en el río Nexapa, encontraron 
que los sitios más contaminados están donde 
se cultivan flores y hortalizas, pues se aplican 
mayores cantidades de agroquímicos y se 
riega con agua residual. 


Las concentraciones de las formas químicas 
del nitrógeno en los ríos varían de acuerdo con 
las condiciones específicas de cada caso. Por 
ejemplo, para el río Texcoco, la concentración 
promedio de nitrato (NO,) fue igual a 0.56 mgL- 
1 (0.13 como N-NO,) (Navarro-Garza y Pérez- 
Olvera, 2005) y de nitrógeno total de 21 a 60 
mgL” (Guzmán et al., 2007). López et al. (2007), 
en un estudio realizado en el cauce principal 
de la cuenca Lerma-Chapala, reportan valores 
de N-NH, y N-NO, de 3.97 y 12.92 mgL”, 
respectivamente. La concentración de N-NO, 
en el Lerma-Chapala se aproxima al reportado 
por Burt et al. (2011), para los ríos de Inglaterra 
y Gales, donde se alcanzan valores de 10 mgL”. 
Concentraciones menores a las del río Texcoco 
fueron reportados por Santiago-Fragoso y 
Sandoval-Manrique (2001) para el río Cuautla, 
quienes encontraron valores de 0.01 a 0.13 
mgL” de N-NH,, de 0.005 a 0.1 mgL” de N,,,, y 
de 0.02 a 0.23 mgL” de NT. 

Las concentraciones de nitrógeno en los ríos 
tienen variaciones temporales. Xia et al. (2002) 
encontraron que el N-NH, y el NIT mostraron 
un acelerado aumento a través de los años 
1980, 1990, 1997 y 1999, atribuyéndolo a los 
incrementos de la descarga de agua residual 
y a la aplicación de fertilizantes nitrogenados. 
Además, las concentraciones de N-NH,, NIT y 
N-NO, son mayores en la época de sequía que 
en la de lluvias. En el río Amarillo (China), la 
concentración de N-NH, pasó de 0.04 en 1980 
a 0.36 mgL” en 1999, y la del NIT pasó de 1.3 
en 1990 a 4.64 mgL* en 1999. Las acciones que 
se tomen para reducir las concentraciones de 
nitrógeno tendrán su efecto a mediano y largo 
plazos. Así lo demuestra el trabajo de Burt et al. 
(2011), quienes señalan que las concentraciones 
de nitrato en algunas corrientes de agua de 
cuencas dominadas por el agua subterránea, 
por ejemplo el río Frome (Reino Unido), 
siguen aumentando, aunque ya se han tomado 
medidas de control de descargas de aguas 
residuales a los ríos. Con base en lo anterior 
y para tener una idea de lo que ha pasado 
respecto a las concentraciones de nitrógeno en 
el río Zahuapan (México), es muy importante 


Muñoz-Nava et al., Relación entre el nitrógeno y las microcuencas del río Zahuapan, México 


describir espacial y temporalmente, de 1997 
a 2007, las concentraciones de las especies 
de nitrógeno (total, orgánico, amoniacal y 
nitratos) de la corriente de este río, y su relación 
con algunos atributos de las microcuencas, 
población y precipitación pluvial. 


Metodología 


La cuenca del río Balsas pertenece a la región 
hidrológica RH18 de la república mexicana 
(Conagua, 2010). La cuenca del río Balsas está 
dividida en tres subregiones hidrológicas: 
a) cuenca del Alto Balsas, b) cuenca del río 
Cutzamala y el Medio Balsas, y c) cuenca del 
Bajo Balsas. La cuenca del Alto Balsas o de 
los ríos Atoyac y Mixteco se subdivide en las 
siguientes subcuencas: Las Cerradas Orientales, 
Alto Atoyac, Bajo Atoyac, Nexapa, Amacuzac, 
Tlapaneco y Mixteco (Toledo, 2003). La red de 
drenaje es compleja, puesto que la roca de la 
Faja Volcánica Transmexicana, y por lo tanto de 
la subcuenca, tiene una composición química 
poco típica de los cinturones volcánicos 
asociados con los límites convergentes de 
placas (Martínez et al., 2007). En la subcuenca 
del Alto Atoyac hay escurrimientos que forman 
dos corrientes principales: el del río Atoyac y 
el del río Zahuapan el segundo es considerado 
afluente del primero. Este estudio se realizó en 
el río Zahuapan, estado de Tlaxcala (figura 1) 
y su área de captación es mencionada en este 
trabajo como subcuenca. La corriente del río 
Zahuapan se forma de los escurrimientos de 
la sierra de Tlaxco, Bloque de Tlaxcala y de 
la Malintzin. En la sierra de Tlaxco nace el río 
Zahuapan y se forma la corriente El Fondón. 
Del Bloque de Tlaxcala se forman las corrientes 
Totolac y Rojano, y de la Malinztin se forman las 
corrientes Amomoloc, Atenco, Grande y Viejo. 
Al margen del río hay tres manantiales en los 
municipios de Tlaxco, Amaxac y Apetatitlán 
(observación directa realizada en campo). 
Se desconoce cuál es la aportación de estos 
manantiales al caudal del río. En la temporada 
de sequía, un porcentaje importante del caudal 
es aportado por la presa de Atlangatepec y por 


las descargas de aguas residuales. Hay siete 
represas que interrumpen la corriente del río 
y está asentado un total de mil localidades 
en la subcuenca (estimación hecha por los 
autores). El río Zahuapan fluye de norte a sur 
y está dividido en dos segmentos por la presa 
de Atlangatepec (figura 1). El segmento que 
descarga a la presa tiene una longitud de 23 km 
(segmento 1); la del que se une al río Atoyac 
es de 75 km (segmento 2). La pendiente media 
del río es de 0.011 (estimación hecha por los 
autores). En la subcuenca del río Zahuapan 
hay un mosaico de varios tipos de suelo. Los 
que más área cubren son los de tipo fluvisol, 
regosol, cambisol, vertisol, feozem, litosol y 
andosol (Werner 1988). La precipitación media 
anual es de aproximadamente 700 mm. En el 
norte de la subcuenca, la actividad económica 
predominante es la agricultura de temporal; 
en el sur hay una combinación de actividades 
agrícola (de riego y de temporal) e industrial. 
La población humana se concentra en el 
centro-sur de la subcuenca, con densidades 
poblacionales en algunos municipios mayores 
a 1 000 hab km” (estimación de los autores). 

Se utilizaron cuatro fuentes de información 
relacionada con la concentración de especies 
de nitrógeno en el río Zahuapan. En total se 
recopilaron 1 782 datos de nitrógeno orgánico 
(Nor) nitrógeno como amoniaco (N-NH,), 
nitrógeno como nitrito (N-NO,), nitrógeno 
como nitrato (N-NO,), nitrógeno total Kjeldahl 
(NTK) y nitrógeno total (NT). ELNTK es la suma 
de N,,, Y N-NH, y el NT es la suma de NTK, 
N-NO, y N-NO,. Los datos de concentración 
de nitrógeno más antiguos son de 1978 (Sedue, 
1985). En ese trabajo se reportan 14 datos de 
concentración de NTK, para siete puntos de 
muestreo distribuidos en todo el río. Otros 175 
datos registrados en 1997, 35 para cada uno de 
las siguientes especies: Na N-NH,, N-NO,, 
NTK y NT fueron publicados en un estudio 
preliminar del río Zahuapan, realizado por la 
Coordinación General de Ecología del Estado 
de Tlaxcala (CGE, 1997). En dicho estudio se 
hizo el muestreo del agua de la corriente del río 
en 13 puntos, de marzo a agosto de 1997. Esos 
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Figura 1. Localización del área de estudio. Las números 1 y 2 indican los segmentos en los que está dividido 


el río Zahuapan por la presa de Atlangatepec. 


13 puntos de muestreo abarcaron un segmento 
de 39 km en la parte sur del río Zahuapan. La 
tercera fuente de información fue la Comisión 
Nacional del Agua, de esta dependencia se ob- 
tuvieron 503 datos de concentración de cuatro 
especies de nitrógeno: 54 de Na 63 de N-NH,, 
142 de N-NO,, 140 de N-NO,, 54 de NTK y 50 
datos de NT, para los años de 1999 a 2007. Estos 
datos corresponden a tres puntos de muestreo 
localizados en el norte, centro y sur del río, que 
se encuentran en la Red Nacional de Monitoreo 
de esta misma dependencia. La cuarta fuente 
de datos fue un proyecto de investigación 
desarrollado por la Universidad Autónoma de 
Tlaxcala, financiado por SIZA Conacyt (Clave 
de Registro TLAX-2004-C01-27). En ese estudio 
se llevaron a cabo muestreos mensuales (por 
los autores) de septiembre de 2006 a agosto de 
2007, en 22 puntos del río Zahuapan. Los puntos 
de muestreo se seleccionaron considerando los 
principales tributarios del río. Se obtuvieron 
en total 1 090 datos de concentración de N,,,, 
N-NH,, N-NO,, NTK y NT, cada una de estas 


especies tuvo 218 datos. Las determinaciones 
analíticas de las formas de nitrógeno fueron 
hechas por instituciones oficiales de acuerdo 
con las normas mexicanas vigentes en esos 
momentos. Aunque se han hecho enmiendas 
a las normas, éstas no han modificado los 
principios químicos de determinación, por 
lo que los datos aquí utilizados pueden 
considerarse comparables entre sí. Más aún, la 
información de la primera, segunda y cuarta 
fuente fue generada por el Laboratorio de lo 
que hoy es la Coordinación General de Ecología 
del Gobierno del Estado de Tlaxcala. 

El cálculo de los atributos: área, pendiente, 
porcentajes del área con vegetación natural 
(comprende varios tipos de bosques, al igual 
que chaparral y pastizales), agricultura y áreas 
erosionadas, así como densidad de población 
en cada microcuenca se realizó de la siguiente 
manera: las capas de las curvas de nivel a 10 
m (INEGL 1999), que cubren el estado de 
Tlaxcala, se procesaron en ArcMap'% versión 
9.3 (ESRI?, 2008) para generar el modelo digital 
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de elevación del tipo TIN (Triangulated Irregular 
Network). El programa TAS (Terrain Analysis 
System) versión 2.0.9 (Lindsay, 2005), al ser 
alimentado con el modelo digital de elevación 
y los puntos de descarga, produjo los grid's, los 
cuales se procesaron en ArcMap'Y, para generar 
los polígonos de las microcuencas (figura 2). 
En ArcMap'Y, los polígonos de los atributos 
publicados en mapas digitales por INEGI (1995) 
para el estado de Tlaxcala fueron cortados 
por los polígonos de las microcuencas, para 
finalmente calcular las áreas con este mismo 
programa. Se calculó la tasa de la densidad de 
población en un nivel de microcuenca, para 
el lapso de 1990 a 2010, mediante regresiones 
lineales. Para ello se utilizaron los censos y 
conteos de población del INEGI de 1995 a 2010 
(INEGI, 2011a-d). 

Las concentraciones de las especies químicas 
de nitrógeno se analizaron en el nivel de punto 
de muestreo y de microcuenca (agrupando los 
puntos de muestreo que fueron incluidos en 
los polígonos de los parteaguas, es decir, que 


pertenecen a una microcuenca). Se graficaron 
las concentraciones de NT, NTK, N-NH,, 
Nos N-NO, y N-NO, contra los años de 1978 
a 2007, para visualizar las tendencias en ese 
periodo en la subcuenca del río Zahuapan. 
Se realizó estadística descriptiva, que incluye 
la media, mediana, moda, mínimo, máximo, 
varianza, desviación estándar, coeficiente de 
variación, varianza y Kurtosis, empleando el 
programa Statistica versión 8.0 (StatSoft, Inc., 
2007). Se determinaron las tendencias de las 
concentraciones de NT, NTK, N-NH, y N,, y en 
el nivel de microcuencas. Con base en que en el 
altiplano mexicano se definen con claridad dos 
épocas del año, una de lluvias y otra de sequía, 
se agruparon los datos de las concentraciones 
de nitrógeno (NT, NTK, N-NH, No. dl N-NO, y 
N-NO,) con este criterio. Se utilizó la definición 
de meses de sequía y de lluvias publicada por 
la Conagua (2011). En los años donde hubo 
datos, la época de sequía se formó con los 
meses de los años subsecuentes (la época de 
sequía en el altiplano es de noviembre a mayo). 
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Figura 2. Microcuencas del río Zahuapan. Los números indican los puntos de muestreo. 
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Se calcularon los coeficientes de correlación 
relacionando las concentraciones de las especies 
químicas de nitrógeno de los años 2006 y 2007 
en siete puntos de muestreo contra los datos 
de precipitación pluvial mensual y diario; para 
este último caso, la registrada un día antes 
de la fecha de muestreo en siete estaciones 
meteorológicas. También se determinaron 
coeficientes de correlación formando pares de 
datos entre las concentraciones de nitrógeno en 
el nivel de punto de muestreo y de microcuenca. 
Para determinar el nivel de significancia de los 
valores de R (coeficiente de correlación), se 
utilizó la prueba de t (McBean y Rovers, 1998). 
En el nivel de microcuenca se realizaron análisis 
de regresión múltiple de las concentraciones de 
nitrógeno de las formas N-NO, = A, + A,*N- 
NA, + AN. y NINO, =4,+ A NT +A, UN- 
NH, + AN ,, donde los coeficientes A, son 
los valores de los coeficientes de regresión 


múltiple lineal. 
Resultados 


El resultado de la ejecución del programa 
TAS fueron los polígonos de las microcuencas 
etiquetadas con lasletras de Aa L (figura2). Estos 
polígonos son los límites de los parteaguas de 
las doce microcuencas. Los puntos de muestreo 
están enumerados del 1 al 27, empezando en 
la parte alta y terminando en la parte baja de 
la subcuenca. Algunos puntos de muestreo, 
principalmente de años anteriores a 2006- 
2007, quedaron dentro de los polígonos de las 
microcuencas H, K y L. El área, la pendiente, 
los porcentajes del área con vegetación natural, 
con agricultura y erosionadas, así como la 
densidad de población de las microcuencas 
del río Zahuapan, se encuentran en el cuadro 
1. Las áreas van de 18.1 a 529.4 km”. La más 
pequeña fue la microcuenca A y la más grande 
la L. El área total de las doce microcuencas 
representó el 48% de la superficie del estado 
de Tlaxcala. La microcuenca A, donde se 
encuentra el nacimiento del río, resultó tener 
la mayor pendiente, con un valor de 11.5%. 
La microcuenca D, que se encuentra después 


de la presa de Atlangatepec, tuvo la menor 
pendiente, con 2.2%. Respecto al porcentaje 
del área cubierta por vegetación natural (la 
vegetación natural comprende varios tipos de 
bosques, así como chaparral y pastizales), la 
microcuenca A tuvo la mayor cobertura, con 
un porcentaje de 92.1%. Las microcuencas 
E y F tuvieron una cubierta vegetal natural 
que no alcanza el 0.01% de su superficie. En 
el nivel de subcuenca, la del río Zahuapan 
tuvo un porcentaje de cubierta vegetal igual 
a 30.5%. En relación con el porcentaje del área 
dedicada a la agricultura, la microcuenca D 
tuvo el porcentaje más alto, 92.2%, seguida de 
H, E J y C. Por subcuenca, el porcentaje del 
área dedicada a la agricultura tuvo un valor de 
61.2%. La microcuenca E sobresalió por haber 
tenido el mayor porcentaje de área erosionada 
con 31.1%; por otro lado, la microcuenca 
A tuvo solamente 0.5%. La densidad de 
población promedio, calculada con los censos 
y conteos de 1990 a 2010, fue de 307.2 hab-km?” 
en la subcuenca. Las microcuencas l, | y H 
tuvieron las mayores densidades de población 
promedio, alcanzando en este lapso valores de 
1 153.8, 964.4 y 734.6 hab-km”, respectivamente. 
Con datos de población del censo de 2010, las 
densidades de estas mismas microcuencas 
alcanzaron valores de 1 422.7, 1 721.7 y 902.13 
hab:km?. Las pendientes de las regresiones 
lineales de las densidades de población contra 
el tiempo indican que estas tres microcuencas 
tuvieron las mayores tasas de incremento de 
densidad de población, con valores de 28.9, 
74.0 y 16.0 hab:-km”*año”, respectivamente, en 
el periodo comprendido de 1995 a 2010. 
Respecto a las concentraciones de las 
especies de nitrógeno en las microcuencas 
(figura 3a), las líneas de las tendencias de 
las concentraciones de NT, NTK, N-NH, y 
N , muestran que la microcuenca A tuvo los 
menores promedios, sigue un incremento en 
la microcuenca B, y una disminución en la C 
y D. A partir de esta última microcuenca hubo 
un aumento sostenido hasta la microcuenca 
L, que es donde se encuentra la salida del río 
de la subcuenca. Las concentraciones en la 
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microcuenca L, comparadas con las de la D, 
se multiplicaron por 7 para NT y NTK, por 
10 para N-NH, y por 8 para N,,. En cuanto a 
N-NO, y N-NO,, las líneas muestran que las 
concentraciones subieron y bajaron a lo largo 
del río, teniendo una alza de 1.1 y 0.13 mgL- 
!, para N-NO, y N-NO,, respectivamente, 
entre estas dos microcuencas (D y L). En 
forma puntual (figura 3b), las tendencias de 
las curvas de concentraciones de NT, NTK 
y N-NH, muestran claramente cinco picos 
en los puntos de muestreo 3, 10, 13, 18 y 25 
(microcuencas B, E, H y J), y dos caídas en los 
puntos 7 y 22-23 (microcuencas E y L). En los 
puntos de muestreo 25 a 27 (no pertenecen a 
ninguna microcuenca considerada en este 
trabajo), antes de la unión del río Zahuapan 
con el río Atoyac, se observaron las mayores 
concentraciones promedio del N,,,. En el punto 
17 se advirtió el pico de mayor concentración 
de N-NO,. 

Las concentraciones promedio de las formas 
de nitrógeno de los datos agrupados con base 
en los años (figura 3c) fueron 10.1, 8.71, 5.87, 
4.23, 0.45 y 0.07 mgL”, para NT, NTK, N-NH, 
Nory N-NO, y N-NO,, respectivamente. Las 
curvas de concentraciones de NT, NTK y 
N-NO, mostraron el mayor valor en el año 2001 
y las de N-NH,, Ne y N-NO, en 2002; a partir 
de esos años, las curvas tendieron a disminuir. 
Estas máximas concentraciones fueron 29.0 
mgL” para NT y NTK, y 6.2 mgL” para N-NO.. 
Las concentraciones de N-NH,, N,,, y N-NO, 
alcanzaron valores máximos iguales a 9.9, 6.2 y 
5.0 mgL”, de manera respectiva. Con excepción 
del N-NH,, las concentraciones de las demás 
especies químicas en 2007 fueron menores, 
comparadas con las de 1997. Los datos de NTK 
de 1978 tuvieron el menor valor de Kurtosis, 
esto es, tuvieron la distribución más platicúrtica; 
los datos de NTK de 2000 tuvieron un valor de 
Kurtosis cercano a cero, es decir, tuvieron una 
distribución mesocúrtica, o más semejante a la 
distribución normal; los datos de N-NO, para 
2000 tuvieron el valor de Kurtosis más alto, lo 
que indica que presentaron la distribución más 
leptocúrtica. En relación con el coeficiente de 


variación, los datos de N-NO, mostraron los 
mayores valores, alcanzando 300% en al año 
2000; por otro lado, las concentraciones de 
N-NH, tuvieron los menores coeficientes de 
variación de alrededor de 30% en el año 2002. 

En el altiplano mexicano se definen clara- 
mente dos épocas del año: una de lluvias y 
otra de sequía. Al agrupar los datos con base 
en el criterio época del año, se encontró que 
las concentraciones de NT, NTK, N-NH, y 
No; fueron mayores en el estiaje que en la 
época de lluvias; el caso contrario ocurrió 
para las concentraciones de N-NO, y N-NO,. 
Las especies que tuvieron las mayores dife- 
rencias de concentración entre la época de 
estiaje y de lluvias fueron N-NO, y N-NO,, que 
alcanzaron 55 y 220%, respectivamente. En el 
periodo que comprende el estudio, se observó 
que las mayores concentraciones de NT, NTK y 
N-NH, ocurrieron en la época de estiaje, entre 
los años 2002-2003 y 2003-2004 (figura 3d). 
Estas mismas especies de nitrógeno tuvieron 
las menores concentraciones en las épocas de 
lluvias de los años 2003 y 2007, pero en este 
último año, las concentraciones se redujeron 
aún más. Respecto al N-NO,, las mayores 
concentraciones se observaron en las épocas 
de lluvias de 1997, 2001, y para el N-NO,, la 
mayor concentración se presentó en la época 
de lluvias de 2002. El promedio mensual de las 
concentraciones entre 1997 y 2007 mostraron 
que las mayores concentraciones promedio de 
NT, NTK y N-NH, se presentaron en diciembre. 
Los promedios más bajos de estas mismas 
especies se tuvieron en agosto. Para el caso del 
N-NO, y N-NO,, las mayores concentraciones 
promedio se registraron en julio. 

De los 60 coeficientes de correlación que 
se realizaron, entre las concentraciones y 
las precipitaciones pluviales mensuales, 29 
resultaron mayores que 0.5 o menores que 
-0.5. Sin embargo, al aplicar la prueba de t, 24 
coeficientes de correlación resultaron con un 
nivel de significancia al 0.05, y siete coeficientes 
con un nivel de significancia de 0.01. En cuanto 
alas correlaciones entre las concentraciones y la 
precipitación pluvial diaria, 14 de 52 resultaron 
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Figura 3. Concentraciones promedio de las especies de nitrógeno en el río Zahuapan: a) por microcuenca, 


b) por puntos de muestreo, c) por año y d) época del año. 


mayores que 0.5 o menores que -0.5. Con la 
prueba de t, sólo ocho coeficientes resultaron 
con un nivel de significancia al 0.05 y cero a un 
nivel de 0.01. 

En el nivel de punto de muestreo, los 
coeficientes de correlación de NT-N-NH, 
y NEN,,, resultaron en 26 y 23 de 27 casos, 
respectivamente, mayores que 0.5. Con la 
prueba de t se obtuvieron 23 y 18 coeficientes, 
con un nivel de significancia de 0.05 y 20, y 12 
al 0.01 para estas mismas especies químicas. 
Para los pares N-NH¿-N,,., N,,¿N-NO, y NT- 
N-NO, se obtuvieron valores de coeficientes 
de correlación mayores que 0.5 o menores que 


0.5, en 9, 7 y 7 casos, respectivamente. Con la 
prueba de t, el número de coeficientes con un 
nivel de significancia al 0.05 fueron 7, 3 y 5; 
para 0.01 se redujeron a 4, 1 y 3 coeficientes, 
respectivamente. Por microcuenca, los resul- 
tados también mostraron coeficientes de corre- 
lación con un nivel de significancia de 0.05 y 
0.01 entre NT-N-NH, y NI-N,,,. El análisis 
de regresión múltiple de la forma N-NO, = 
A, + A¿N-NH, + AN e mostró coeficientes 
de correlación mayores que 0.5 en tres de las 
12 microcuencas; con la forma N-NO,= A, + 
ANT 4 A NENE + AN se obtuvieron 
coeficientes de correlación mayores que 0.5 en 
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nueve de las 12 microcuencas. Los valores de 
los coeficientes de las ecuaciones de regresión 
y los valores de R se encuentran en el cuadro 2. 


Discusión 


Los promedios de concentración de NT en 
las doce microcuencas no fueron mayores 
a 15 mgL”, este valor es el límite máximo 
promedio diario establecido en la NOM-001- 
ECOL-1996 (Semarnat, 2003), para protección 
de vida acuática en ríos. La microcuenca L es 
la más próxima a este límite, con un valor de 
11.5 mgL”. 

Al comparar cada dato con este límite, se 
encuentra que 38 de ellos fueron mayores que 
éste. En esta microcuenca se sobrepasó en más 
ocasiones (ocho veces) este límite y como se 
encuentra a la salida de la subcuenca del río, esto 
hace suponer que el NT se va acumulando de 
la cabecera a la salida de la subcuenca. También 
estos promedios de NT fueron menores a los 
reportados por Guzmán et al. (2007), quienes 
encontraron valores que oscilaron entre 21 y 
60 mgL” para el río Texcoco. En los criterios 
ecológicos CE-CCA-001/89 (Sedue, 1989), se 
establece como nivel máximo para nitrógeno 
amoniacal (N-NH,) una concentración de 0.06 


mgL" para protección de vida acuática de agua 
dulce. La concentración promedio de N-NH, 
en la subcuenca fue de 4.54 mgL”. De los 315 
datos, 310 fueron mayores que este límite. 
En la microcuenca C se registraron los cinco 
datos menores a este valor. La concentración 
promedio de N-NO, igual a 0.32 mgL” en la 
subcuenca del río Zahuapan fue mayor a lo 
reportado por Navarro-Garza y Pérez-Olvera 
(2005) para el río Texcoco, quienes reportaron 
concentraciones menores de 0.13 mgL”, pero 
fueron más bajas que lo reportado para el 
río Frome, que alcanzó valores de 10 mgL”, 
debido a que este río tiene entradas de agua 
subterránea (Burt et al., 2011). 

En la figura 3a se ve con claridad un 
descenso de las concentraciones de las 
especies de nitrógeno, aquí consideradas, 
de la microcuenca B a la D. Entre estas 
dos microcuencas se encuentra la presa de 
Atlangatepec, en este embalse hay plantas 
acuáticas; así como a lo largo de unos 2 km 
en el cauce del río después de la cortina de la 
presa. Este tramo del río funciona como un 
pequeño humedal. Para el caso particular de 
estas dos microcuencas, el NT disminuyó de 
8.3 a 4.3 mgL” (48%), el NTK de 8.1 a 4.1 mgL” 
(49%), el N-NH, de 4.8 a 2.2 mgL” (54%), el 


Cuadro 2. Coeficientes de las ecuaciones de regresión múltiple entre las concentraciones de las especies de nitrógeno por 


microcuenca. R es el coeficiente de correlación. 


Ecuación de regresión Ecuación de regresión 
i N-NO,= A, +A,N-NH, + 4,N,,, N-NO,=A,+A,NT+A,'N-NH,+A/N,,, |  R 
Microcuenca R 
A, A, A, A, A, A, A, 
A 0.08 0.07 0.05 0.34 0.14 0.43 -0.44 -0.39 0.67 
B 0.30 -0.01 -0.01 0.25 0.29 0.05 -0.05 -0.05 0.32 
E 0.20 -0.01 0.06 0.39 0.13 0.27 -0.28 -0.22 0.62 
D 0.35 0.03 -0.10 0.57 0.34 0.04 -0.01 -0,14 0.60 
E 0.18 0.07 0.02 0.57 0,12 0.28 -0.23 -0.27 0.73 
F 0.18 -0.02 0.14 0.38 0.13 0.25 -0.26 -0.14 0.59 
G 0.27 -0.03 0.09 0.45 0.18 0.25 -0.27 -0.18 0.63 
H 0.03 0.01 0.10 0.31 0.01 0.21 -0.20 -0.13 0.54 
I -0.21 -0.06 0.37 0.60 =0:12 0.75 075 -0.64 0.92 
J 0.47 -0.08 0.11 0.38 0.06 0.71 -0.73 -0.67 0.87 
K 0.19 -0.02 0.05 0.38 0.16 0.11 -0.13 -0.07 0.49 
L 0.34 -0.01 -0.01 0.23 0.32 0.06 -0.06 -0.06 0.33 
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N,;, de 3.3 a 1.9 mgL” (42%) y el N-NO, de 
0.24 a 0.23 mgL” (4%). No fue posible estimar 
la baja del nitrito por carecer de datos. Estos 
porcentajes de diminución, con excepción de 
los nitratos, se aproximan al 50% de remoción 
de N reportado por Kunz et al. (2011b), para 
la presa Itezhi-Tezhi, localizada en el bajo río 
Kafue, el mayor tributario del río Zambezi 
(Zambia), y por Kunz et al. (2011la), para 
el lago Kariba (Zambia-Zimbawe), ambos 
grandes cuerpos de agua localizados en clima 
tropical. La desnitrificación y acumulación 
en sedimentos manera 
importante en la remoción de nitrógeno, 
pudiendo alcanzar 70% (Kunz et al., 2011a). 
En los embalses, las especies químicas del 


contribuyen de 


nitrógeno están sujetas a reducción del 
nitrato a nitrógeno gaseoso en los sedimentos, 
absorción del amonio por las plantas acuáticas 
y fijación del amonio en las arcillas (Stumm y 
Morgan, 1996). Estos procesos en la presa de 
Atlagantepec posiblemente sean los causantes 
de la disminución de las concentraciones de 
las especies químicas del nitrógeno. 

Conlos datos del cuadro 1 hubo correlación 
positiva con un nivel de significancia al 0.01 
entre el área acumulada de las microcuencas 
en dirección de la corriente del río y las 
concentraciones de las especies de nitrógeno 
para NT y N-NH,. La correlación positiva indi- 
ca que el área acumulada y las concentracio- 
nes de estas especies de nitrógeno tuvieron la 
misma tendencia creciente de la cabecera a la 
salida de la subcuenca. Sin embargo, para el 
caso del N-NO,, la correlación resultó igual 
que -0.415, pero no fue significativa. Respecto 
a la pendiente del terreno de las microcuencas 
y las concentraciones de las especies químicas 
medidas en la corriente del río, se encontró 
que no hay correlación. El mayor valor de 
correlación se obtuvo para la pendiente contra 
el N-NO,, con una r = -0.404, la cual no es 
significativa. De la misma forma, los datos de 
las concentraciones de nitrógeno no tuvieron 
relación con las áreas de cubierta vegetal, 
agricultura y erosión en las microcuencas. 
Este resultado difiere de lo reportado por 


otros autores. Por ejemplo, Tesoreiro et al. 
(2009) señalan una relación entre el nitrato y 
la agricultura, Burt et al. (2011) y Alexander et 
al. (2008) mencionan la relación existente en su 
estudio entre las concentraciones de nitrógeno 
y la agricultura. Por otro lado, los datos de las 
densidades de población en las microcuencas 
y las concentraciones de NT, N-NH, y N-NO, 
mostraron correlaciones positivas a un nivel 
de significancia de 0.05 con la prueba de t. Con 
este análisis se deduce que las concentraciones 
de NT y N-NH, están correlacionadas con 
el área acumulada de las microcuencas 
conforme a la dirección de la corriente y la 
densidad de población, está última también 
tuvo correlación con el N-NO.. 

El comportamiento de las concentraciones 
de las especies químicas del nitrógeno en los 
años considerados mostró que de 1997 a 2001 
hubo un incremento, y de este año a 2007, 
una disminución de NT, NTK y N-NO),; esto 
mismo ocurrió para N-NH,, N,,, y N-NO,, 
pero el cambio de tendencia ocurre en 2002. 
Xia et al. (2002) reportaron para el río Amarillo 
un incremento de 0.04, en 1980, a 0.36 mgL” 
de N-NH, en 1999, es decir, una tasa anual 
de 0.0168; mientras que en el río Zahuapan, 
el incremento para N-NH, es de 3.86 en 1997 
a 8.11 mgL* en 2002, lo cual arrojó una tasa 
anual de 0.85. De 2002 a 2005, la tasa anual 
de disminución de la concentración de 
esta especie química fue de 0.798, con una 
diferencia entre 2007 y 1997 de 0.26 mgL-. 
Como el área, la pendiente del terreno de 
las microcuencas, la cubierta vegetal, la 
agricultura y áreas erosionadas, no han 
cambiado en las microcuencas en este lapso 
y la población ha mantenido un aumento 
constante, la única explicación del cambio de 
tendencia de las concentraciones son factores 
externos que promueven la salida de nitró- 
geno del río. Una posible salida de nitrógeno 
del río puedeser las plantas de tratamiento, que 
tratan el agua residual de algunos municipios 
antes de descargar al río Zahuapan, y por 
condiciones ambientales de la corriente del 
río, que promueven desnitrificación, es decir, 
la transformación de NO”, = N,. 
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También Xia et al. (2002) reportaron para 
el río Amarillo que las concentraciones de 
N-NH, NIT y N-NO, son mayores en la época 
de sequía que en la de lluvias. Este mismo 
resultado se obtuvo para NT, N-NH, y N,,, en 
el río Zahuapan. Sin embargo, el N-NO, y el 
N-NO, tuvieron las mayores concentraciones 
en la época de lluvias, en los años 1997 y 
2001 para el primero, y 2002 para el segundo. 
En específico, las mayores concentraciones 
de N-NO, y N-NO, se presentaron en julio, 
es decir, en plena época de lluvias, sin que 
esto signifique que exista correlación entre la 
precipitación pluvial y las concentraciones 
de nitrógeno, debido a que no se obtuvieron 
coeficientes de correlación significativos. En 
las microcuencas I y J se notó claramente la 
diferencia de las concentraciones entre estas 
dos épocas del año. No se vio de manera 
contundente qué variable considerada en este 
trabajo tuvo el mayor grado de influencia 
en este resultado. Sin embargo, en esas dos 
microcuencas, las densidades de población 
fueron las más altas comparadas con las otras. 

Por último, con base en Snoeyink y Jenkins 
(1990), y Stumm y Morgan (1996), se puede 
plantear que la cadena de transformaciones 
que tiene el nitrógeno en los puntos mues- 
treados considerados en este trabajo del 
río Zahuapan es el siguiente: N,, (del agua 
residual) > NH; > NO; > NO; => N,. Esto 
es, que no hay otras entradas de nitrato al 
río, como por ejemplo debidas a las aguas 
subterráneas y que la salida es sólo por N,, 
del cual no se tienen mediciones, por lo que 
no es considerado en el análisis de regresión 
múltiple. Esta cadena de trasformaciones del 
nitrógeno se fundamenta con el resultado del 
análisis de regresión múltiple, que considera 
al NO, como una variable dependiente del 
NT, Mos y NH,, con base en que se obtuvieron 
coeficientes de correlación mayores que 0.5 en 
las microcuencas A, C, D, E, E G, H, I y J. En 
estas dos últimas microcuencas se obtuvieron 
valores de coeficientes de correlación iguales 
a 0.92 y 0.87, y las mayores densidades de 
población (1 153 y 964 hab-km7). 


Conclusiones 


Con la realización de este trabajo, empleando 
atributos de las microcuencas del río Zahuapan 
y datos históricos de concentraciones de NT, 
Nory N-NH,, N-NO, y N-NO, en la corriente 
del río, se concluye que: 


e Las concentraciones promedio de NT 
no excedieron, en ninguna de las doce 
microcuencas, el límite máximo promedio 
diario (15 mgL”) establecido en la NOM- 
001-ECOL-1996 (Semarnat, 2003), para 
protección de vida acuática en ríos. 

e En toda la subcuenca del río Zahuapan, 
la concentración de N-NH, fue mayor a 
0.06 mgL”, que es el límite establecido en 
los criterios ecológicos CE-CCA-001/89 
(Sedue, 1989) para la protección de vida 
acuática de agua dulce. Sin embargo, las 
concentraciones promedio de NT fueron 
menores y las de N-NO, mayores a las 
del río Texcoco. Las concentraciones de 
nitrógeno como nitrato en el río Zahuapan 
fueron menores a las del río Frome. 

e Es posible que las disminuciones de 
concentración de las especies de nitrógeno 
en las microcuencas B, C y D hayan sido un 
efecto de la presa de Atlangatepec y de la 
vegetación acuática que se encuentra en el 
cauce del río Zahuapan. 

+ Se encontró correlación entre las concen- 
traciones de NT y N-NH, con el área acu- 
mulada de las microcuencas, conforme 
a la dirección de la corriente del río y 
la densidad de población, esta última 
también tuvo correlación con el N-NO,. Las 
variables: pendiente del terreno, cubierta 
vegetal, agricultura y áreas erosionadas de 
las microcuencas, no tuvieron correlación 
con ninguna de las especies de nitrógeno 
aquí consideradas. 

e La tendencia creciente de las concen- 
traciones del nitrógeno cambió de una 
tendencia en aumento a una de descenso 
en los años 2001 y 2002. Para el caso de 
N-NH, la tasa de aumento anual, calculada 
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de 1997 a 2002, en el río Zahuapan fue de 
0.85. Pero también hubo una tasa anual 
de disminución de esta especie, de 2002 a 
2005, igual a 0.8. 

+ Para el NT, N-NH, y N,,,, se encontró que 
sus concentraciones disminuyeron en la 
época de lluvias, pero no así para las de 
N-NO, y N-NO,, que aumentaron en esa 
misma época del año. 

e El análisis de regresión múltiple, consi- 
derando al NO, como una variable 
dependiente del NT, Mo y NH,, mostró 
altos valores de coeficientes de correlación. 
Esto sugiere que el NO, en el río Zahuapan 
no depende de otras fuentes, sólo de las 
aguas de origen residual. 
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Abstract 


MUÑOZ-NAVA, H., SUÁREZ-SÁNCHEZ, J., VERA-REYES, A., BAUMANN, ]J., 
CARREÓN-COCA, M.F. € MONTEALEGRE-MUÑOZ, R.G. Relationship between 
nitrogen and the Zahuapan River microbasins, Mexico. Water Technology and Sciences (in 
Spanish). Vol. IV, No. 4, September-October, 2013, pp. 47-61. 


Agricultural and wastewater discharges are significant sources of nitrogen in streams and 
rivers. The aim of this study was to spatially and temporally describe the concentrations of 
nitrogen species in the Zahuapan River and their relationships with some of the attributes of 
the microbasins, including area, slope, vegetation cover, agriculture, soil erosion, population 
and pluvial precipitation. The polygons for the Zahuapan River microbasins were defined 
using the TAS program. The results showed that nitrogen concentrations diminished in the 
microbasins located near the Atlangatepec dam. The concentrations were highest during the 
years 2001 and 2002. Concentrations of NT, NTK, N-NH,, and Norg increased during the dry 
months while N-NO3 and N-NO2 decreased during this period of the year. The relationship 
between concentrations and monthly precipitation resulted in a greater number of significant 
coefficient correlations than for daily precipitation. The correlation coefficients for nitrogen 
species were significant at the 0.05 level. The multiple regression analysis N-NO, = A, + 
ANT + A,'N-NH, + AN, showed coefficients greater than 0.5 in 9 to 12 microbasins. 
This study concluded that the Atlangatepec dam plays a role in the concentration of nitrogen 
in the Zahuapan River. 


Keywords: nitrogen, Zahuapan River, microbasins. 
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